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We resort to the concept of turns to provide a geometrical representation of the action of any 
lossless multilayer, which can be considered as the analogous in the unit disk to the sliding vectors in 
Euclidean geometry. This construction clearly shows the peculiar effects arising in the composition 
of multilayers. A simple optical experiment revealing the appearance of the Wigner angle is analyzed 
in this framework. 
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PACS numbers: 230.4170 Multilayers, 120.5700 Reflection, 120.7000 Transmission, 000.3860 Mathematical 
methods in physics 



I. INTRODUCTION 

The search for mathematical entities that could de- 
scribe physical phenomena has always been a top priority. 
For example, the need to describe a direction in space, 
combined with the use of geometry to approach physical 
problems, brought forth the concept of a vector. 

The idea that the complex numbers have a geometrical 
interpretation as vectors lying in a plane, led to Hamilton 
to introduce quaternions with the aim of being useful for 
the analysis of three-dimensional space . The price to 
be paid is that the composition of quaternions is not com- 
mutative. Soon after that, it became clear that rotations 
can be advantageously represented by unit quaternions. 

A notion closely related to Hamilton treatment is that 
of turns 2]. The turn associated with a rotation of axis 
n and angle $ is a directed arc of length t?/2 on the great 
circle orthogonal to n on the unit sphere. By means of 
these objects, the composition of rotations is described 
through a parallelogramlike law: if these turns are trans- 
lated on the great circles, until the head of the arc of 
the first rotation coincides with the tail of the arc of the 
second one, then the turn between the free tail and the 
head is associated with the resultant rotation. Hamilton 
turns are thus analogous for spherical geometry to the 
sliding vectors in Euclidean geometry. 

In recent years many concepts of geometrical nature 
have been introduced to gain further insights into the 
behavior of layered media. The algebraic basis for these 
developments is the fact that the transfer matrix associ- 
ated with a lossless multilayer is an element of the group 
SU(1, 1), which is locally isomorphic to the Lorentz group 
SO(2, 1) in (2+1) dimensions. This leads to a natural and 
complete identification between reflection and transmis- 
sion coefficients and the parameters of the corresponding 
Lorentz transformation 0, Q . 

In an appealing paper, Juarez and Santander || devel- 
oped a generalization of Hamilton turns to the Lorentz 
group, while Simon, Mukunda, and Sudarshan 0, [?J 
worked out an equivalent algebraic approach for SU(1, 1), 



in which they introduce a noncommutative geometrical 
addition for these hyperbolic turns that reproduces the 
composition law of the group; i. e., for both the reflec- 
tion and transmission coefficients. The latter coefficients 
seem to be (almost) ignored in the literature. The goal of 
this paper is precisely to show how this formalism affords 
a very intuitive image of multilayer optics. 



Of course, since the amplitude and phase of a light 
beam can be conveniently represented by a (Euclidean) 
vector, the properties of layered media have been rep- 
resented by graphical constructions. These visualization 
tools are commonly in use in thin-film design and include 
the Smith chart, the admittance diagram 0, the (re- 
flectance) circle diagrams , and the vector method ^3 > 
among others. However, the formalism of turns is intrin- 
sic and does not rely on any specific representation for the 
light waves. In other words, once the turn for a multilayer 
is known, his action on any light state is fully determined. 
We emphasize that, although, these techniques can be 
used for quantitative calculations, they cannot compete 
nowadays with modern software design, and their great 
value is in the visualization of the characteristics of a 
multilayer. 



This paper is organized as follows: In Section 2 we 
present some details of how the action of any multilayer 
can be seen as a geometrical motion in the unit disk. Ev- 
ery one of these motions can be decomposed in terms of 
two reflections, which justifies the idea of turn introduced 
in Section 3, where their composition law is also intro- 
duced via a parallelogram law, in close analogy with what 
happens for sliding vectors in Euclidean geometry. The 
noncommutative character of this law leads to interest- 
ing phenomena, such as the appearance of extra phases 
in the composition of multilayers, which is examined in 
Section 4, giving a simple though nontrivial example that 
illustrates clearly how this geometrical scenario works in 
practice. 
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FIG. 1: Wave vectors of the input [E a and E s ] and output 
[Ei~^ and £^ + '] fields in a multilayer sandwiched between two 
identical semi-infinite ambient and substrate media. 



II. MULTILAYER ACTION IN THE UNIT DISK 

The theory of reflection and transmission of light by 
stratified planar structures is of wide interest in optics. 
A considerable amount of theoretical work has been done 
on this topic, a detailed discussion of which can be found 
in a number of books 

BE IH IH 13 Nowadays, 
the standard method uses a matrix representation of the 
field in each medium, as pioneered by Abeles In 
this paper, we follow the elegant approach developed by 
Hayficld and White 

We deal with a stratified structure that consists of a 
stack of plane-parallel layers sandwiched between two 
semi-infinite ambient (a) and substrate (s) media that 
we shall assume to be identical, in order to simplify as 
much as possible the calculations. Hereafter all the me- 
dia are supposed to be lossless, linear, homogeneous, and 
isotropic. We choose the Z axis perpendicular to the 
boundaries and directed as in Fig. 1. 

A monochromatic linearly polarized plane wave falls 
from the ambient making an angle 9q with the normal 

to the first interface and with an amplitude E a + ^ . The 
electric field is either in the plane of incidence (p po- 



larization) or perpendicular to the plane of incidence (s 
polarization). Since the multilayer has two input chan- 
nels, we consider as well another plane wave of the same 
frequency and polarization, and with amplitude E^ , 
incident from the substrate at the same angle 6q. This 
guarantees that all the multiple reflected and transmitted 
waves superpose with the same wave vectors. The over- 
all output fields in the ambient and the substrate will be 
denoted E a ~^ and e[ + \ respectively. In the common ex- 
perimental situation, there is no light incident from the 
substrate [Es = 0], although there are relevant situa- 
tions in which both input fields are present (e. g., at the 
beam splitter in a Michelson interferometer). 

The field amplitudes at each side of the multilayer are 
related by the linear relation 



E { - ] 



M„ 



(+) 

s 

(-) 



(2.1) 



where the multilayer transfer matrix M as can be shown 
to be 



1/T as Rasl^as 
Ras/Tas V^as 





a 








a* 



(2.2) 



Here the complex numbers R as and T asi which can be 
expressed as 

R as = \R as \ exp(ip), T as = \T as \ exp(ir), (2.3) 

are, respectively, the overall reflection and transmission 
coefficients for a wave incident from the ambient. Be- 
cause |i? QS | 2 + |T as | 2 = 1, we have det M as = \a\ 2 -\(3\ 2 = 
1, and then M as belongs to the group SU(1, 1). 

In Ref. we have proposed viewing the multilayer 
action in a relativisticlike framework, giving a formal 
equivalence between the fields in Eq. (|2.1|) and space- 
time coordinates in a (2 + l)-dimensional space. These 
coordinates verify (x ) 2 ~ (x 1 ) 2 — (x 2 ) 2 — 1, which de- 
fines a unit two-sheeted hyperboloid characteristic of the 
group SO(2, 1). If one uses stereographic projection tak- 
ing the south pole as projection center, the upper sheet of 
the unit hyperboloid is projected into the unit disk, the 
lower sheet into the external region, while the infinity 
goes to the boundary of the unit disk. 

The geodesies in the hyperboloid are intersections with 
the hyperboloid of planes passing through the origin. 
Consequently, hyperbolic lines are obtained from these 
by stereographic projection and they correspond to cir- 
cle arcs that orthogonally cut the boundary of the unit 
disk. 

In many instances (e.g., in polarization optics |13|| ') we 
are interested in the transformation properties of field 
quotients rather than the fields themselves. Therefore, it 
seems natural to consider the complex numbers 
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(2.4) 



The action of the multilayer given in Eq. (|2.2() can be then 
seen as a function z a — f(z s ) that can be appropriately 
called the multilayer transfer function |l8j . 

From a geometrical viewpoint, this function defines a 
transformation of the complex plane, mapping the point 
z s into the point z a , according to 



z a = $[M as ,z s ] = 



(3* + a*z s 
a + (3z s 



(2.5) 



which is a bilinear or Mobius transformation. One can 
check that the unit disk, the external region and the 
unit circle remain invariant under the multilayer action. 
Note also that, when no light impinges from the sub- 
strate z s = and then z a = R as - It is worth mentioning 
that this approach is quite general, since it provides the 
transformed of any point in the unit disk for every values 
of the input fields E ( a +) and E ( s ~ ] . 

To classify the multilayer action it proves convenient to 
work out the fixed points of the mapping; that is, the field 
configurations such that z a = z s = Zf in Eq. (|2.5|l |2Cj : 



z f 



*[M r . 



whose solutions are 



Zf = 2/3 ^ Im(a) ± ^ Tr ( M -)] 2 - 4 ] 



(2.6) 



(2.7) 



When [Tr(M as )] 2 < 4 the multilayer action is elliptic and 
it has only one fixed point inside the unit disk. Since in 
the Euclidean geometry a rotation is characterized for 
having only one invariant point, this multilayer action 
can be appropriately called a hyperbolic rotation. 

When [Tr(M as )] 2 > 4 the action is hyperbolic and it 
has two fixed points, both on the boundary of the unit 
disk. The geodesic line joining these two fixed points re- 
mains invariant and thus, by analogy with the Euclidean 
case, this action will be called a hyperbolic translation. 

Finally, when [Tr(M as )] 2 = 4 the multilayer action is 
parabolic and it has only one (double) fixed point on the 
boundary of the unit disk. 

Here we will be concerned only with the case 
[Tr(M a;! )] 2 > 4, since it is known that any element of 
SU(1, 1) can be written (in many ways) as the product 
of two hyperbolic translations 7] . The axis of the hyper- 
bolic translation is the geodesic line joining the two fixed 
points. A point on the axis will be translated to another 
point, a (hyperbolic) distance |2l| 



C = 2 In 



l - {Tr(M as ) + v/[Tr(M as )] 2 - 4} 



(2.8) 



along the axis. 



III. HYPERBOLIC TURNS AND THEIR 
COMPOSITION 

In Euclidean geometry, a translation of magnitude £ 
along a line 7 can be seen as the product of two reflec- 
tions in any two straight lines orthogonal to 7, separated 




FIG. 2: Representation of the sliding turn T-j / 2 in terms of 
two reflections in two lines Fi and F2 orthogonal to the axis 
of the translation 7, which has two fixed points Zf\ and 2/2. 
The transformation of a typical off axis point z 3 is also shown. 



a distance £/2. This idea can be translated much in the 
same way to the unit disk, once the concepts of line and 
distance are understood in the hyperbolic sense. In con- 
sequence, any pair of points z\ and zi on the axis of the 
translation 7 at a distance £/2 can be chosen as intersec- 
tions of Ti and T2 (orthogonal lines to 7) with 7. It is 
then natural to associate to the translation an oriented 
segment of length (/2 on 7, but otherwise free to slide 
on 7 (see Fig. 2). This is analogous to Hamilton's turns, 
and will be called a hyperbolic turn T 7 ^/ 2 - 

Note that using this construction, an off-axis point 
such as z s will be mapped by these two reflections 
(through an intermediate point z- lnt ) to another point 
z a along a curve equidistant to the axis. These other 
curves, unlike the axis of translation, are not hyperbolic 
lines. The essential point is that once the turn is known, 
the transformation of every point in the unit disk is au- 
tomatically established. 

Alternatively, we can formulate the concept of turn 
as follows. Let M as be a hyperbolic translation with 
Tr(M as ) positive (equivalently, Re(a) > 1). Then, M as 
is positive definite and one can ensure that its square root 
exists and reads as 



1 



y/2[Bc(a) + 1] 



a + 1 (3 
13* a* + 1 



(3.1) 



This matrix has the same fixed points as M as , but the 
translated distance is just half the induced by M as ; i.e., 



C(M as ) = 2C( V / M^T). 



(3.2) 



This suggests that the matrix \/M as can be appropri- 
ately associated to the turn T 7) £ / 2 that represents the 
translation induced by M as . Therefore, we symbolically 



FIG. 3: Composition of two hyperbolic turns T 



T 



71X1/2 



and 



72X2/2 



by using a parallelogramlike law when the axes 71 



and 72 of the translations intersect. 



write 



FIG. 4: Composition of two multilayers represented by Her- 
mitian matrices Hi and H2. H2 maps the point z 3 — — R2 into 
the origin, while Hi maps the origin into z a = R\. We show 
also the associated turns Ti and T2, as well as the resulting 
one T( 12 ) obtained via the parallelogram law. The composite 
multilayer H1H2 transforms the point z s into z a . The data of 
the corresponding multilayers are given in the text. 



T. 



7.C/2 



(3.3) 



One may be tempted to extend the Euclidean com- 
position of concurrent vectors to the problem of hyper- 
bolic turns. Indeed, this can be done quite straightfor- 
wardly |5(. Let us consider the case of the composition 
of two of these multilayers represented by the matrices 
Mi and M2 (for simplicity, we shall henceforth omit the 
subscript as) of parameters Ci and (2 along intersecting 
axes 71 and 72, respectively. Take the associated turns 
T 7li ^/2 and T 72j( ; 2 /2 and slide them along 71 and 72 until 
they are "head to tail" . Afterwards, the turn determined 
by the free tail and head is the turn associated to the re- 
sultant, which represents thus a translation of parameter 
£ along the line 7. This construction is shown in Fig. 3, 
where the noncommutative character is also evident. 

In Euclidean geometry, the resultant of this parallel- 
ogram law can be quantitatively determined by a direct 
application of the cosine theorem. For any hyperbolic 
triangle with sides of lengths £1 and C2 that make an an- 
gle 6, we take the expression from any standard book on 
hyperbolic geometry pll ] 



cosh Q = cosh Ci cosh (2 + sinh £1 sinh £ 2 cos 8, (3.4) 

where 9 is the angle between both sides. Moreover, for 
future use we quote that the (hyperbolic) area VL of the 
geodesic triangle is 



tan(fi/2) = 



tanh(Ci/2) tanh(£ 2 /2) sin6> 
1 - tanh(Ci/2) tanh(C 2 /2) cos( 



(3.5) 



IV. APPLICATION: REVEALING THE 
WIGNER ANGLE IN THE UNIT DISK 

To show how this formalism can account for the exis- 
tence of peculiar effects in the composition of multilayers, 
we address here the question of the Wigner angle in the 
unit disk and propose a simple optical experiment to de- 
termine this angle. 

The Wigner angle emerges in the study of the com- 
position of two noncollinear pure boosts in special rel- 
ativity: the combination of two such successive boosts 
cannot result in a pure boost, but renders an additional 
pure rotation, usually known as the Wigner rotation [2^] 
(sometimes the name of Thomas rotation [23[ 01 IS & ^ so 
used). In other words, boosts are not a group. 

To fix the physical background, consider three frames 
of reference K, K' and K" . Frames K-K' and K'-K" 
have parallel respective axes. Frame K" moves with uni- 
form velocity V2 with respect to K' , which in turn moves 
with velocity Vi relative to K . The Lorentz transforma- 
tion that connects K with K" is given by the product 
Li(vi)l-2(v2), which can be decomposed as 



Li( Vl )L 2 (v 2 ) = L (12) (v)R(tf), 



(4.1) 



where one must be careful in operating just in the same 
order as written in Eq. (|4.1I) . In words, this means that 
an observer in K sees the axes of K" rotated relative to 
the observer's own axes by a Wigner rotation described 
by R(^). The explicit expression for the axis and angle 
of this rotation can be found, e.g. in Ref. [z| and will be 
worked out below from the perspective of the multilayer 
action in the unit disk. 
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First, we observe that any matrix M S SU(1, 1) can 
be expressed in a unique way in the form 



M = HU, 



(4.2) 



where H is positive definite Hermitian and U is unitary. 
One can check by simple inspection that the explicit form 
of this (polar) decomposition reads as [25| 



M 



HU 



1/\T\ R*/\T\ 
R/\T\ 1/\T\ 

exp(— it) 
exp(ir 



(4.3) 



The component H is equivalent to a pure boost, while U 
is equivalent to a spatial rotation. 

It is clear from Eq. Ij4.3[l that [Tr(H)] 2 > 4, so it repre- 
sents a hyperbolic translation. Moreover, one can check 
that its associated fixed points are diametrically located 
on the unit circle and so, the axis of this translation is 
precisely the diameter joining them. By writing now 

R = tanh(C/2) exp(ip), T = scch(C/2) exp(ir), 

(4.4) 

one can easily check that the matrix H in Eq. I|4.3|l trans- 
forms the origin into the complex point R, that is, 



$[H,0] = i?, 



$[i-r\i?i = o. 



(4.5) 



In complete analogy with Eq. (j4.1|l we compose now 
two multilayers represented by Hermitian matrices Hi 
and H2 (that is, with zero transmission phase lag t\ = 
72 = 0) and we get, after simple calculations [2fiT | 



HiH 2 



H(i2)U 



1/IT, 



(12)1 



R 



R(12)/\T(12)\ 

exp(-i*/2) 

exp(i*/2) 



(12)/ \ T (12)\ 
V|7(12)| 



(4.6) 



where 



^(12) - 



R\ + i?2 

1 + R{R 2 



T (12) 



\TiT 2 \ 



1 + R\R 2 ' 



(4.7) 



- = arg[T ( i 2) ] 



= arg(l + R X R* 2 ), 



and the subscripts 1 and 2 refer to the corresponding 
multilayer. The appearance of an extra unitary matrix in 
Eq. (|4.6[) is the signature of a Wigner rotation in the mul- 
tilayer composition and, accordingly, the Wigner angle 
is just twice the phase of the transmission coefficient of 
the compound multilayer. 

To view this Wigner angle in the unit disk, let z s be 
the point in the substrate that is transformed by the 
multilayer H2 into the origin, and let z a be the result of 
transforming the origin by Hi. According to Eq. I|4.5|l . 
one has 



-R2 



0. 



$[Hi,0] = Rl 



(4.8) 
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FIG. 5: The same as in Fig. 4 but now the resulting turn T( 12 ) 
has been slided to three different positions along the axis. The 
corresponding points are also transformed by H1H2. All the 
geodesic triangles plotted have the same hyperbolic area *F 



Consider now the (geodesic) triangle defined by the 
points z s , O and z a in Fig. 4. The general formula 1)3. 5|) 
gives for this triangle 



(4.9) 



which confirms the geometric nature of this Wigner an- 
gle, since it can be understood in terms of the area (or 
equivalently, the anholonomy) of a closed circuit j23> I23 . 

According to the ideas developed in Section 3, we can 
reduce the multilayers Hi and H2 to the associated turns, 
represented by arrows in Fig. 4. The "head to tail" 
rule applied to Ti and T2 (for simplicity, we omit in 
the subscripts of these turns the corresponding param- 
eters) immediately gives the resulting turn T( 12 ) . How- 
ever, note that, if we follow the formal prescription shown 
in Eq. (|3.3|) and ascribe Ti 1— * and T2 1— > VH2, we 
conclude that the composition law imposes 



T 



(12) 



VH1H2 



(4.10) 



All these results are independent on the position of 
the turn. In fact, in Fig. 5 we have put the turn T(i2) 
in different positions along the axis 7. In every position, 
we have drawn two radii passing through the head and 
the tail of T( i2 ) and taken on them twice the hyperbolic 
distance from the origin. The pairs of points obtained in 
this way (such as z' s and z' a ) are transformed precisely by 
H1H2. In other words, H1H2 can be decomposed in many 
ways as the composition of two Hermitian matrices and 
every geodesic triangle z' s O z' a has the same hyperbolic 
area 'F 

It seems pertinent to conclude by showing an exper- 
imental implementation of the data shown in Fig. 4. 
To this end, we first recall 0] that for a single plate 
of refractive index nj and thickness dj imbedded in air 
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H, 



FIG. 6: Scheme of two Hermitian multilayers Hi and H2. The 
compound multilayer H1H2 obtained putting together these 
two components induces a Wigner rotation of angle 2r. 



(no = 1) illuminated with a monochromatic light of wave- 
length in vacuo A, a standard calculation gives the trans- 
fer matrix M(yo with the following reflection and trans- 
mission coefficients: 



OjO 



r 0j [l - exp(-i25j)} 



T, 



().;() 



1 - rL exp(- 



tojtjo exp(- 



428 j) ' 



■iSi) 



(4.11) 



1 - rjy cxp(-i2Sj) ' 



where roj and toj are the Fresnel reflection and trans- 
mission coefficients at the interface Oj (which applies to 
both p and s polarizations by the simple attachment of 
a subscript p or s) and 6j is the plate phase thickness 



Sj = ——rijdjCOsOj, 
A 



(4.12) 



9j being the angle of refraction in the layer. The transfer 
matrix for the coherent addition of m of these layers is 



on a glass substrate (with refractive index 713 = 1.5 and 
thickness 1^3 = 1.3 mm), and imbedded in air. The light 
has a wavelength in vacuo of Ao = 546 nm and falls from 
the ambient at normal incidence. Such a simple system 
could be manufactured with standard evaporation tech- 
niques. 

We have performed a computer simulation of the per- 
formance of this multilayer Hi using a standard package, 
obtaining T x = 0.9055 and R ± = 0.3736 - 0.2014i, which 
in turn gives Ty = and p\ = —0.4944 rad. 

Our second multilayer H2 is a symmetric system 
formed by two films of zinc sulphide and thickness d± = 
40 nm separated by a spacer of air with a phase thickness 
S 5 = 3.707 rad. For this subsystem we have T 2 = 0.9399 
and i?2 = 0.3413z, and therefore T2 = and p 2 = tt/2 
rad. 

When these two multilayer are put together by the 
marked points in Fig. 6, the resulting one has a trans- 
mission phase lag of r = —0.1361 rad, which is just 
half the area of the geodesic triangle z s Oz a in Fig. 4, as 
predicted by the theory. 

In summary, we expect that the geometrical approach 
presented here will be an interesting tool for representing 
in a graphical way the multilayer action. Moreover, the 
composition law of these turns allows for a clear under- 
standing of the nontrivial effects appearing in the com- 
position of multilayers. 

We stress that the benefit of this approach lies not in 
any inherent advantage in terms of efficiency in solving 
problems in layered structures. Rather, we expect that 
turns could provide a general and unifying tool to ana- 
lyze multilayer performance in an elegant way that, in 
addition, is closely related to other fields of physics. 



M = ]J M 



(4.13) 



As shown in Fig. 6, we take as the first multilayer Hi 
the lossless system formed by two thin films, one of zinc 
sulphide (with refractive index ri\ = 2.3 and thickness 
d\ = 80 nm) and the other of cryolite (with refractive 
index 122 = 1.35 and thickness di — 104 nm), deposited 
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